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The excitation of two particle-two hole final states in neutrino-nucleus scattering has been ad¬ 
vocated by many authors as the source of the excess cross section observed by the MiniBooNE 
Collaboration in the quasi elastic sector. We analyse the mechanisms leading to the appearance of 
these final states, and illustrate their significance through the results of accurate calculations of the 
nuclear electromagnetic response in the transverse channel. A novel approach, allowing for a consis¬ 
tent treatment of the amplitudes involving one- and two-nucleon currents in the kinematical region 
in which the non relativistic approximation breaks down is outlined, and its preliminary results are 
reported. 
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I. INTRODUCTION 

Experimental studies of neutrino-nucleus interactions 
carried out over the past decade um have provided am¬ 
ple evidence of the inadequacy of the Relativistic Fermi 
Gas Model (RFGM), routinely employed in event gen¬ 
erators, to account for both the complexity of nuclear 
dynamics and the variety of reaction mechanisms—other 
than single nucleon knock out—contributing to the ob¬ 
served cross section. 

A striking manifestation of the above problem is the 
large discrepancy between the predictions of Monte Carlo 
simulations and the double differential charged current 
(CC) quasi elastic (QE) cross section measured by the 
MiniBooNE Collaboration using a carbon target [3]. 

As pointed out by the authors of Ref. j5], improving 
the treatment of nuclear effects, which turns out to be 
one of the main sources of systematic uncertainty in the 
oscillation analysis [6], will require the development of 
a comprehensive and consistent description of neutrino- 
nucleus interactions, validated through extensive compar¬ 
ison to the large body of electron-nucleus scattering data 

[BIB]. 

The main difficulty involved in the generalisation of the 
approaches successfully employed to study electron scat¬ 
tering to the case of neutrino interactions stems from the 
fact that, while the energy of the electron beam is fixed, 
in neutrino scattering the measured cross section results 
from the average over different beam energies, broadly 
distributed according to a flux $. Therefore, a measure¬ 
ment of the energy of the outgoing charged lepton in a 
CC QE interaction does not specify the energy transfer to 
the nuclear target, which largely determines the reaction 
mechanism. As shown in Refs. p[ [10], the MiniBooNE 
double differential cross section corresponding to a spe¬ 
cific muon energy bin turns out to receive comparable 
contributions from different mechanisms, which must be 
all taken into account in a consistent fashion. 

Many authors have suggested that the excess CC QE 


cross section observed by the MiniBooNE collaboration 
is to be ascribed to the occurrence of events with two 
particle-two hole final states, not taken into account by 
the RFGM employed for data analysis [3 HUH2. The 
description of these processes within a realistic model 
of nuclear dynamics requires that all mechanisms lead¬ 
ing to their occurrence—Initial State Correlations (ISC) 
among nucleons in the target nucleus, Final State Corre¬ 
lations (FSC) between the struck nucleon and the spec¬ 
tator particles, and interactions involving two-nucleon 
meson-exchange currents (MEC)- be included. Within 
the Independent Particle Model (IPM) of the nucleus, 
however, correlations are not taken into account, and two 
particle-two hole final states can only be excited through 
the action of two-body operators, such as those involved 
in the definition of MEC. 

In this paper, we analyse the mechanisms leading to 
the appearance of two particle-two hole final states in the 
response of interacting many-body systems, and argue 
that the interference between amplitudes involving one- 
and two-nucleon currents plays an important role. This 
feature clearly emerges from the results of a calculation of 
the transverse electromagnetic response of 4 He and of the 
corresponding sum rule of 12 C, computed using state-of- 
the-art models of the nuclear hamiltonian and currents, 
within the Green’s Function Monte Carlo (GFMC) com¬ 
putational scheme [13] . 

In view of the extension of our study to the kinematical 
regime in which the non relativistic approximation is no 
longer applicable, we also outline a novel approach, based 
on a generalisation of the factorisation ansatz , underly¬ 
ing the spectral function formalism. This scheme, allow¬ 
ing for a consistent treatment of one- and two-nucleon 
current contributions, appears to be quite promising for 
applications to neutrino scattering. 

The structure of the nuclear cross section, as well as its 
expression in terms of longitudinal and transverse struc¬ 
ture functions, are reviewed in Section [n] while Sec¬ 
tion m describes the theoretical approaches, based on 
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nuclear many-body theory, developed to study electron- 
nucleus scattering. In Section III A we discuss the non 
relativistic regime and the results of GFMC calculations, 
while in Sections |III B| and |III C| we derive the explicit 
expression of the two particle-two hole contribution to 
the cross section obtained from our approach, and dis¬ 
cuss the preliminary results of its application. Finally, 
in Section |IV| we summarise our findings and state the 
conclusions. 


II. NUCLEAR CROSS SECTION AND 
RESPONSE FUNCTIONS 


In the one-photon-exchange approximation, the double 
differential electron-nucleus cross section can be written 
in the form 


d 2 o 

dE e 'dtt e > 


or E e 
E t 




(1) 


where k e = (E e , k e ) and k e i = (E e >, k e /) are the four- 
momenta of the incoming and outgoing electrons, respec¬ 
tively, a = 1/137 is the fine structure constant, dCl e > is 
the differential solid angle in the direction specified by 
k e /, and q = k e — k e ' = (w,q) is the four momentum 
transfer. 

The lepton tensor L pi/ is completely determined by 
lepton kinematics, while the nuclear response is described 
by the tensor IF//'/ defined as 


where 


Al = 



At 


1 q 2 2 

+ tan 

2 q 2 


Oe 
2 ’ 


( 6 ) 


and (do/dVl^M = [acos(9 e /2)/AE e sin 2 (0 e /2)] 2 is the 
Mott cross section. 

The L and T structure functions can be readily ex¬ 
pressed in terms of the components of the response ten¬ 
sor of Eq. ([2]). Choosing the 2 -axis along the direction 
of the momentum transfer one finds 


Rl = Wf (7) 

RT=t (<% - q -^f) W% = W? + Wf . (8) 

ij =1 q 

Note that the above expressions are completely 
general, and describe processes involving both one- and 
two-nucleon current operators. 

It follows from Eqs. ([l]) and ([ 2 ]) that the nuclear cross 
section and response functions can be written as a sum 
of contributions corresponding to different hadronic final 
states |IV). Consider, for example, the case of QE scat¬ 
tering, in which the final state particles are nucleons only. 
For a carbon target we find 

\N) = \ n B : p) , | U C» ,\ w B lP n) ,\ 10 Be,pp)..., (9) 


WT(q.w) -£>l JH\N)(N\JZ\0)6W(P o + q-p N ), 

N 

( 2 ) 


where |0) and |iV) denote the initial and final hadronic 
states, the four-momenta of which are Pq = (Eq, po) and 
P N = (EpjjPjy). The nuclear current can be written as 
a sum of one- and two-nucleon contributions, according 
to (see, e.g., Ref. mi) 

= (3) 

i j>i 


The current j? describes interactions involving a single 
nucleon. In the QE sector, it can be expressed in terms of 
the measured vector form factors [15] . The two-nucleon 
contribution ;)/•, on the other hand, accounts for pro¬ 
cesses in which the beam particle couples to the currents 
arising from meson exchange between two interacting nu¬ 
cleons. 

Equation 0 can be rewritten in terms of two response 
functions, denoted i?z,(q, w) and ^^(q, w), describing in¬ 
teractions with longitudinally (L) and transversely (T) 
polarised photons, respectively. The resulting expression 
reads 


d 2 a 

dE' e d£l e 



A L (|q|,w,6» e )R L (|q|,o;) 


+ A T (\q\,u,8 e )R T (\q\,w) , 


( 4 ) 

( 5 ) 


where the residual nucleus can be in any bound state. 

The states | N) are usually classified according to the 
number of nucleons excited to the continuum, and re¬ 
ferred to as one particle-one hole (lplh), two particle- 
two hole (2p2h), etc. In Eq. ([ 9 ]), \ n B,p) and | n C,ri) 
are lplh states, while | w B 1 pn) and \ 10 Be,pp) are 2p2h 
states. 

Neglecting the contributions of final states involving 
more than two nucleons in the continuum, the cross sec¬ 
tion can be written as 

do = doi p i h + do 2p2 h oc L pu {W^ lh + W^ 2h ) . (10) 

We recall that, in scattering processes involving inter¬ 
acting many-body systems, 2p2h final states can be pro¬ 
duced through the action of both one- and two-nucleon 
currents]/] However, in order for the matrix element of a 
one-body operator between the target ground state and a 
2p2h final state to be non vanishing, the effects of dynam¬ 
ical nucleon-nucleon (NN) correlations must be included 
in the description of the nuclear wave functions. 

Correlations give rise to virtual scattering between nu¬ 
cleons in the target nucleus, leading to the excitation 
of the participating particles to continuum states. The 


1 It should be kept in mind that lplh final states can also be 
excited by both one- and two-nucleon currents. 
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ISC contribution to the 2p2h amplitude arises from pro¬ 
cesses in which the beam particle couples to one of these 
high-momentum nucleons. The FSC contribution, on the 
other hand, originates from scattering processes involv¬ 
ing the struck nucleon and one of the spectator particles, 
that also result in the appearance of 2p2h final states. 


III. MANY-BODY THEORY OF THE NUCLEAR 
RESPONSE 

As discussed in the previous Section, the calculation 
of the nuclear response requires the evaluation of the 
transition amplitudes (0| J^\N), involving both one- and 
two-nucleon current operators, as well as all possible fi¬ 
nal states. The initial state can be accurately described 
within the framework of non relativistic many-body the¬ 
ory using realistic models of the nuclear hamiltonian, 
strongly constrained by nucleon-nucleon scattering data 
and nuclear phenomenology. The final state and the cur¬ 
rent operator, on the other hand, depend on momentum 
transfer, and their calculation in the kinematical region 
in which the non relativistic picture breaks down neces¬ 
sarily implies additional assumptions. 


A. Non relativistic regime 

The approach based on the GFMC computational 
scheme provides a suitable framework to carry out ac¬ 
curate calculations of a variety of nuclear properties in 
the non relativistic regime, typically corresponding to 
|q| < 500 MeV (for a recent review of Quantum Monte 
Carlo methods for nuclear physics see, e.g., Ref. USD- 

Valuable information on the L and T responses can be 
obtained from their Laplace transforms, also referred to 
as Euclidean responses, defined as 
/»00 

Et,l{ q,r)= / due~ ulT RT,L( q,w) ■ (11) 

Juiel 

The lower integration limit w e 1 = q 2 /2M a, Ma being 
the mass of the target nucleus, is the threshold of elastic 
scattering—corresponding to the | N) = |0) term in the 
sum of Eq. § -the contribution of which is excluded. 

Within GFMC, the Euclidean responses are evaluated 
from 


^(q,r) = (0|p*(q)e-^-^)V(q)|0} 

-|(0|p(q)|0)| 2 e-^, (12) 


and 

= (0|j^(q)e" (H ” Bo)T j T (q)|0) 
-|(0|j T (q)|0)| 2 e-^, (13) 

where p(q) and jT(q) denote non relativistic reductions 
of the nuclear charge and transverse current operators, 
respectively El- 


Note that, although the states | N) ^ |0) do not ap¬ 
pear explicitly in Eqs. (121 and (13), the Euclidean re¬ 


sponses include the effects of final state interactions (FSI) 
of the particles involved in the electromagnetic interac¬ 
tion, both among themselves and with the spectator nu¬ 
cleons. 

The Euclidean responses at r = 0 are directly related 
to the sum rules of the L and T responses, obtained from 
w-integration after removing the trivial energy and mo¬ 
mentum dependence associated with the nucleon form 
factor 1171 : 


C T l f°° 

S'T.idql) = , nP in2 x 12 / dujR T ,L( q,w) . 

NbWqeJJ Ju., 


(14) 


In the above equation, G^(Qq E ) is the electric proton 
form factor evaluated in quasi elastic kinematics, i.e. at 
Qqe = q 2 -w^ EJ with W QE = (a/ q 2 +to 2 -to), where m 
is the proton mass. The coefficients appearing in Eq. (14) 
are defined as 


C L = 


1 

Z 


C T = 


z^i 


Nut 


(15) 


where Z is the proton charge, N = A — Z is the number 
of neutrons and /i p and fi n are the proton and neutron 
magnetic moments, respectively. 

The inversion of the Laplace transform, needed to re¬ 
trieve the energy dependence of the responses, is long 
known to involve severe difficulties. A groundbreaking re¬ 
sult has been recently reported by the authors of Ref. [TB] , 
who exploited the maximum entropy technique to obtain 
the L and T responses of 4 He. 

Figure [I] shows the breakdown of the transverse re¬ 
sponse of H4e at |q| = 500 MeV—computed within the 
approach of Ref. m — into one-nucleon current, two- 
nucleon current and interference contributions. Note that 
the quantity displayed in the figure is normalised dividing 
by the squared proton form factor. 

It clearly appears that including the two-nucleon cur¬ 
rents leads to a sizeable enhancement of the response, 
and that the large positive contribution of the interfer¬ 
ence term peaks at energy loss to < wqe- This feature 
is a direct consequence of nucleon-nucleon correlations, 
neglected in the mean field approach. The agreement be¬ 
tween the GFMC results and the data of Ref. [T3| turns 
out to be remarkably good. 

The extension of the procedure employed to obtain the 
4 He response to heavier nuclei, such as carbon, is still 
out of reach of the available computational capabilities. 
However, valuable information can be extracted from the 
analyses of the sum rules. 

The results of numerical calculations of the carbon 
St (|q|), displayed in Fig. i clearly show that interfer¬ 
ence terms provide a sizeable fraction of the sum rule. 
At momentum transfer |q| > 300 MeV, their contribu¬ 
tion turns out to be comparable to—in fact even larger 
than—that obtained squaring the matrix element of the 
two-nucleon current. 
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FIG. 1. (color online) Transverse response function of 4 He, 
obtained within the approach of Ref. [IS!. The shaded area 
shows the results of the full calculation, with the associated 
uncertainty arising from the inversion of the Euclidean re¬ 
sponse, while the dotted line has been obtained including the 
one-nucleon current only. The dot-dash line represents the 
response computed neglecting the interference term, the con¬ 
tribution of which is displayed by the dashes. The data are 
taken from Ref. m- 



q [MeV] 

FIG. 2. (color online) Sum rule of the electromagnetic re¬ 
sponse of carbon in the transverse channel. The dashed line 
shows the results obtained including the one-nucleon current 
only, while the solid line corresponds to the full calculation. 
The dot-dash line represents the sum rule computed neglect¬ 
ing the interference term, the contribution of which is dis¬ 
played by the dotted line. The results are normalised so that 
the dashed line approaches unity as |q| — > oo. Monte Carlo 
errors bars are not visible on the scale of the figure. 


B. Relativistic regime: the factorisation ansatz 

The results of Figs. |T| and [2] clearly point to the need for 
a consistent treatment of correlations and MEC, within a 
formalism suitable for application in the kinematical re¬ 
gion in which the non relativistic approximation is known 
to fail. This section describes the derivation of the ap¬ 
proach based on factorisation of the nuclear matrix el¬ 
ements. For ease of presentation, we will consider the 


response of uniform and isospin symmetric nuclear mat¬ 
ter. However, the generalisation to atomic nuclei does 
not involve any substantial problems. 

The effects of ISC on the nuclear cross section at large 
momentum transfer can be taken into account using the 
spectral function formalism [19] [20]. The conceptual 
framework underlying this approach is provided by the 
impulse approximation (IA), i.e. the assumption that at 
momentum transfer such that |q| _1 -C d, d being the av¬ 
erage separation distance between nucleons in the target 
nucleus, the nuclear cross section reduces to the incoher¬ 
ent sum of cross sections describing scattering processes 
involving individual nucleons. As a consequence, the con¬ 
tribution of the two-nucleon current can be disregarded, 
and the final state | N) of Eq. © can be written in the 
factorized form 


I N) = |p) <g> \n A -i,Pn) ■ (16) 

In the above equation, |p) is the state of a non inter¬ 
acting nucleon carrying momentum p, while |rM.-i,Pn) 
describes the (A — l)-particle spectator system in the 
state n, with momentum p„. Note that, owing to NN 
correlations, |riA-i,Pn) is not restricted to be a bound 
state [see Eq. ©]. 

Within the IA, the contribution to the nuclear cross 
section arising from interactions involving the one- 
nucleon current can be written in terms of the cross 
sections of elementary scattering processes off individual 
nucleons, the momentum (k) and removal energy (E) of 
which are distributed according to the spectral function 
P(k,E) [IS], defined as 

P{k,E) =Y / \(nA-i,Pn\aw\0)\ 2 6(E + E o -E n ) . (17) 

n 

In the above equation, E n is the energy of the (A — 1)- 
nucleon state, and the operator Ok removes a nucleon of 
momentum k from the nuclear ground state. 

The resulting expression of the cross section is [8j 

d<JiA = f d 3 k dE P(k,E) doi . (18) 


Note that P(k,E) describes an intrinsic property of the 
target nucleus, independent of momentum transfer, and 
as such can be safely obtained from non relativistic many- 
body theory. On the other hand, the matrix elements of 
the nucleon current entering the definition of der; can be 
computed using its fully relativistic form. 

Exploiting the Kallen-Lehman representation of the 
two-point Green’s function, the spectral function ap¬ 
pearing in Eq. (18), can be conveniently split into two 
parts, displaying distinctly different energy dependences 
m- The single particle part Pih(k,E), obtained from 
Eq. ( |17[ ) including bound lh states only, exhibits a pole 
at E = — efc, e*, being the energy of a nucleon in the hole 
state of momentum k. The continuum part, on the other 
hand, is smooth, and extends to large values of energy 
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and momentum. Its leading term, corresponding to 2hlp nucleon is excited to a state outside the Fermi sea, can 
states of the residual ( A— l)-particle system in which one be written in the form 


P 2h i P (k,E) = Jd 3 hd 3 h'd 3 p'\$Z hy l 2 0(k F - |h| )0(k F - |h'|)0(|p'| - k F )6(E + e h 


ep' ) 


(19) 


where the integration includes a sum over the indices 
associated with discrete degrees of freedom, and 


®k h ' P ' = <°l{|k) ® Ihh'p 7 )} . 


( 20 ) 


Note that momentum conservation requires that the ex¬ 
pression of h p involve a 5(h 4- h' — p' k). 

As pointed out above, in the presence of ground state 
correlations both parts of the spectral function pro¬ 
vide non vanishing contributions to the cross section of 
Eq. ph. 
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FIG. 3. (color online) Cross section of the process e + 12 C —► 
e' + X at beam energy E e = 961 MeV and electron scattering 
angle 9 e = 37.5 deg, computed using Eq. (181 with the spec¬ 
tral function of Ref. [20]. The solid line shows the results of 
the full calculation, while the breakdown into lplh and 2p2h 
contributions is illustrated by the dot-dash and dashed lines, 
respectively. 

Figure [3] shows the lplh and 2p2h components of the 
electron-carbon cross section arising from ISC. The cal¬ 
culations have been performed at E e = 961 MeV and 
9 e = 37.5 deg, using Eq. (181 with the spectral func¬ 
tion of Ref. 20] and the parametrisation of the nucleon 
form factors of Ref. [22]. The solid line corresponds to 
the results of the full calculation, while the dot-dash and 


dashed lines have been obtained using the pole and con¬ 
tinuum parts of the spectral function, which amounts to 
taking into account only lplh or 2p2h final states, re¬ 
spectively. The distinct energy dependence of the 2p2h 


r 


contribution, providing ~ 10% of the total QE cross sec¬ 
tion, is clearly visible. 

The importance of relativistic effects can be gauged 
comparing the solid and dashed lines of Fig. [4] repre¬ 


senting the carbon cross sections obtained from Eq. (18) 


using relativistic and non relativistic kinematics, respec¬ 
tively. It clearly appears that in a kinematical setup cor¬ 
responding to |q| ~ 585 MeV at w = wqe relativistic 
kinematics sizeably affects both position and width of 
the quasi elastic peak. 
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FIG. 4. (color online) Electron-carbon cross section obtained 
from Eq. ( |18| ) using relativistic (solid line) and non relativistic 
(dashed line) kinematics. The experimental data are from 
Ref. 1251 . 


The factorisation ansatz of Eq. (16) can be readily 


extended to allow for a consistent treatment of the am¬ 
plitudes involving one- and two-nucleon currents. The 
resulting expression is 


\ N ) = Ipp') ® \m A -2, Pm) 


( 21 ) 


where the states |pp') and |toa-2,Pto) describe two non 
interacting nucleons of momenta p and p' and the (A—2)- 
particle residual system, respectively. 

Using Eq. (21), the nuclear matrix element of the two- 


nucleon current can be written in terms of two-body ma¬ 
trix elements according to 

mm = [ d 3 fcd 3 /c / M ro (k,k')(p P , ij';.|kk'), (22) 
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with M m { k, k') given by 


distribution, defined as 


= {(m (4 _ 2) ,p m | ® (kk'|} |0) . (23) 

From the above equations it follows that the evaluation 
of the nuclear transition matrix element involving the 
two-nucleon current reduces to the calculations of the 
nuclear amplitude M m ( k, k') and of the matrix element 
of the current operator between free nucleon states. The 
former, being independent of momentum transfer, can be 
carried out using the non relativistic formalism, while the 
latter does not involve any approximations. 

The connection with the spectral function formalism 
becomes apparent noting that the two-nucleon spectral 
function P(k,k',£i), yielding the probability of remov¬ 
ing two nucleons of momenta k and k' from the nuclear 
ground state leaving the residual system with excitation 
energy E, is defined as [compare to Eq. @] m 


n(Q) = 4 tt|Q| 2 J d 3 K n + Q, ^ (25) 

where K = k + k', Q = (k — k')/2, and 

n(k,k') = J dE P(k, k', E) , (26) 

is shown by the solid line of Fig. [5j Comparison with the 
prediction of the Fermi Gas (FG) model, represented by 
the dashed line, indicates that correlations give rise to a 
sizeable quenching of the peak of the distribution, along 
with the appearance of a high momentum tail. 

C. lplh and 2p2h contributions to the transition 
matrix element 


P(k,k , ,P) =J2\Mm(KV)\ 2 6(E + E 0 -E m ) , (24) 


with M m ( k, k') given by Eq. (23). 


The two-nucleon spectral function of uniform and 
isospin symmetric nuclear matter at equilibrium density 
has been calculated by the authors of Ref. j23] using a 
realistic hamiltonian. The resulting relative momentum 


The extended factorisation ansatz discussed in the pre¬ 
vious section provides a scheme allowing for a clear cut 
identification of the lplh and 2p2h contributions to the 
nuclear cross section in the presence of two-nucleon cur¬ 
rents. 

Let us consider the contribution of 2p2h final states to 
the response tensor of Eq. ([2]) 


Wg 2h = J d 3 h d 3 ti d 3 p d 3 p9{k F - |h| )8(k F - |h'|)0(|p| - k F )9{ |p'| - k F ) 

x ( 0 |J M |hh'pp')(hh'pp'|P'| 0 ) ( i(w + P 0 -P^p/^q + h + lT - p p') , 


(27) 


where q is the total momentum transfer, while h, h' and 
p, p' are the momenta of the hole and particle states, 
respectively. The structure of the current operator, in¬ 
volving one- and two-body terms, can be best understood 
from its momentum space expression 

J^(k!,k 2 ) = J d 3 X 1 d 3 X 2 J^( Xl ,X 2 ) e-dki-r+k^) 

= jf( k i)£(k 2 ) +^(k 2 )<5(ki) +if 2 (k 1 ,k 2 ) , 

(28) 

showing how the total momentum transfer, 
q = ki + k 2 , is shared between the two nucle¬ 
ons involved in the electromagnetic interaction, labeled 
by the indices 1 and 2. 

Within the factorisation scheme, the matrix element 
of the one-nucleon current operator, can be readily eval¬ 
uated inserting a complete set of states describing a non 
interacting nucleon. The resulting expression is 

(0|jf Ihh'pp') = J d 3 k p <k|jf |p>, (29) 



FIG. 5. (color online) Relative momentum distribution of a 
nucleon pair in isospin symmetric nuclear matter at equilib¬ 
rium density. 


with h p defined by Eq. (20). 


The calculation of the matrix element of the two- 
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nucleon current exploits the fact that, in analogy with 
P(k, E), the two-nucleon spectral function of Eq. ( |24| ) 
can be separated into two parts, characterised by their 
analytical structure. The component corresponding to 
bound 2h states of the (A — 2)-nucleon system exhibits a 
pole located at E = — (e*, + e*/), whereas the continuum 
states, the dominant of which is the 3hlp state, give rise 
to a smooth background. 

It follows that, within the factorisation scheme, the 
contribution to WJ/^h arising from amplitudes involving 
only the two-nucleon current is obtained from the 2 h com¬ 
ponent of P(k.k', E), which can be written in the form 

m 

P 2h (k,k',E) = Id 3 hd 3 h'\^:\ 2 6(E + e h + e h/ ), (30) 
x 0(k F - \h\)0(k F - |h'|) . 

In the above equation, <!>££, is related to the overlap be¬ 
tween the target ground state and the 2h state of the 


(A — 2)-nucleon system through 

^; = (o|{|kk')0|hh')}. (si) 

The diagrammatic analysis of the cluster expansion of 
in uniform and isospin symmetric nuclear matter, 
carried out by the authors of Ref. [3 ]Q, shows that only 
unlinked graphs (i.e., graphs in which the points reached 
by the k, k' lines are not connected to one other by any 
dynamical or statistical correlation lines) survive in the 
A —> oo limit, the contributions of linked diagrams being 
of order 1/A. It follows that 

^; = ^;< 5 (h-k)«j(h'-k'), (32) 

where q is the the Fourier transform of the overlap be¬ 
tween the ground state and the 1 h (A — l)-nucleon state, 
the calculation of which is discussed in Ref. m- 

Collecting the above results, we can write the expres¬ 
sion of the response tensor obtained from the extended 
factorisation ansatz as a sum of three contributions. The 
terms involving the squared amplitudes of the matrix el¬ 
ements involving one- and two-nucleon currents can be 
written in terms of the appropriate contributions to the 
one- and two-nucleon spectral functions, according to 


W 2 p 2 M i = Jd 3 kJ dEP 2hlp (k, £)(k|jf |k + q)(k + q|jj'|k) < y(w -E- e, fc+ ,|)0(|k + q| - k F ) (33) 

and 


Kp2h‘22 = j d 3 kd 3 k'd 3 pd 3 p' J dEP 2h (k,k',E)(kk l \j/ 2 \ P p')(pp , \f 12 \kk l ) (34) 

x <5(k + k' + q - p - p')<5(u> - E - e p - e p /)0(|p| - £;f)0(|p'| - k F ) . 


The interference term, on the other hand, involves a 


product of the nuclear amplitudes entering the definition 
of the spectral functions. The resulting expression is 


W^ 2p 2h,i2 = I d 3 k d 3 i d 3 i' d 3 h d 3 h'd 3 p d 3 p'8( h - £)5(h 7 - £') 


<f>f V <k|/ 


rip; 


*f p (k\f 


2 IP . 


(35) 


x (P,p'Ui 2 l^.O (5(h + h'-f q - p - p')5(u + e h + e h ' - e p - e p /)0(|p| - k F )6{ |p'| - k F ) + h.c. . 


Extensive numerical calculations of the electron- 
carbon cross section based on the formalism described in 
this paper are under way. However, they involve a num¬ 
ber of non trivial novel developments, e.g. the deriva¬ 
tion of the two-hole contributions to the nuclear spectral 
function, the discussion of which will require a separate 
publication. 


Figure [6] shows the transverse electromagnetic response 
of carbon at |q| = 570 MeV computed using the carbon 
spectral function of Ref. [20] and approximating the two- 
hole spectral function of carbon with that of uniform nu¬ 
clear matter, at density corresponding to Fermi momen¬ 
tum k F = 221 MeV. Note that this is not the same as 
working within the Fermi gas model. We use overlaps— 








the functions </>£ defined by Eqs. (31) and (32)—obtained 
within the ab initio approach of Refs. HH12QE4], based 
on a realistic nuclear hamiltonian including two- and 
three- nucleon interactions. 

Owing to short range correlations, which move 
strength from the lplh to the 2p2h sector, the result¬ 
ing occupation of the momentum eigenstates is reduced 
by ~ 20%. 

Interactions effects also affect the initial state energies 
of the knocked-out nucleons [19, 20], thus shifting the 
threshold of the two-nucleon current contributions with 
respect to the predictions of the Frmi gas model EZIHH] 

It has to be pointed out that the correlation contribu¬ 
tion to the carbon spectral function of Ref. [30] is ob¬ 
tained from nuclear matter results. Therefore, the use 
of nuclear matter overlaps in the matrix elements of the 
two-nucleon current entering the interference terms ap¬ 
pears to be consistent. 

We have used the fully relativistic expression of the 
two-nucleon current described in Refs. 123 Hi, with the 
same form factors and A-width. 

The solid line of Fig. [6] represents the results of the 
full calculation, whereas the dashed line has been ob¬ 
tained including only the amplitudes involving the one- 
body current. The contributions arising from the two- 
nucleon current are illustrated by the dash-dot and dot¬ 
ted lines, corresponding to the pure two-body current 
transition probability and the interference term, respec¬ 
tively. The latter turns out to be sizable, its contribution 
being comparable to the total two-body current response 
for w < 350 MeV. Although our results still need to be 
improved, and do not include the corrections taking into 
account the effects of FSI, in Fig. [6]we have also included, 
for comparison, the data resulting from the analysis of 
Ref. (231. 


IV. CONCLUSIONS 

We have analysed the mechanisms—correlations in the 
initial and final states and coupling to meson exchange 
currents—leading to the excitation of 2p2h final states in 
the nuclear response to electromagnetic interactions. In 
the non relativistic regime, in which highly accurate cal¬ 
culations consistently taking into account all these mech¬ 
anisms are feasible, our results confirm the findings of 
Ref. [133. 

In the transverse channel, the contribution of processes 
involving the two-nucleon current is sizeable, and extends 
well into the kinematical region corresponding to energy 
transfer to ~ wqe, in which single nucleon knock out is 
dominant. 

The important role played by interference between the 
amplitudes involving one- and two- body currents clearly 
implies that correlation effects must be included in any 
model aimed at describing the nuclear cross section in 
the 2p2h sector. This point was clearly stated, over three 
decades ago, in the pioneering work of Ref. [23], the au- 



FIG. 6. (color online) Electromagnetic response of carbon 
in the transverse channel, at momentum transfer |q| = 570 
MeV. The solid line represents the results of the full calcu¬ 
lation, whereas the dashed line has been obtained including 
only amplitudes involving the one-body current. The contri¬ 
butions arising from the two-nucleon current are illustrated 
by the dot-dash and dotted lines, corresponding to the pure 
two-body current transition probability and the interference 
term, respectively. The experimental data are taken from 
Ref. [26]. 


thors of which also remarked the inadequacy of the treat¬ 
ment of correlations based on lowest-order perturbative 
pion exchange. However, combining a realistic and con¬ 
sistent description of correlations and meson-exchange 
currents in the kinematical region in which the non rel¬ 
ativistic approximation is no longer applicable involves 
serious difficulties. 

To overcome this problem, we have developed a novel 
approach based on the factorisation ansatz underlying 
the spectral function formalism, widely and successfully 
employed to describe the nuclear response in the lplh 
sector. The preliminary results obtained within this ap¬ 
proach, shown in Fig. [3J provide a fairly good descrip¬ 
tion of the measured transverse response of carbon at 
|q| = 570 MeV. 

A comparison between the results of Fig. [6] and the 
GFMC results of Fig. [ljshows distinctive discrepancies in 
both magnitude and energy dependence of the two-body 
current contributions. While part of the disagreement 
is likely to originate from differences in the two-nucleon 
currents employed in Ref. [18], as well as from the non 
relativistic nature of the GFMC calculations, the large 
interference contribution in the region of the quasi elas¬ 
tic peak observed in Fig. [T] may arise from interference 
between amplitudes involving the one- and two-body cur¬ 
rents and lplh final states. A careful analysis of these 
terms, that were found to be sizable in the pioneering 
work of Ref. [22], is being carried out, and will be dis¬ 
cussed elsewhere. 

The main assumption implied in the factorisation of 
the 2p2h final states is the treatment of the knocked out 
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nucleons as free particles, which amounts to neglecting 
their interactions, both among themselves and with the 
spectator nucleons. Antisymmetrization under exchange 
between any of the outgoing particles and the spectators 
is also disregarded. 

The factorized nuclear transition amplitudes involv¬ 
ing the one-nucleon current can be corrected—to include 
the effects of final state interactions in the quasi elas¬ 
tic sector—using an extension of the spectral function 
formalism, as discussed in Ref. [50]. The resulting modi¬ 
fications lead to (i) a shift in energy transfer of the differ¬ 
ential cross section, arising from interactions between the 
knocked out nucleon and the mean field of the recoiling 
nucleus, and (ii) a redistribution of the strength from the 
quasi free peak to the tails, resulting from rescattering 
processes. Theoretical studies of electron-nucleus scat¬ 
tering suggest that in the kinematical region relevant to 
the MiniBooNE analysis the former mechanism—which 
does not involve the appearance of 2p2h final states— 
provides the dominant contribution, and can be taken 
into account through an optical potential [3lj . 

The corrections to the factorized amplitudes involving 
the two-nucleon current also include interactions between 
the two knocked out particles. A careful analysis of these 
processes is certainly needed. However, the results of 
Shen et al, who carried out an accurate calculation of the 
neutrino-deuteron cross section over a broad kinematical 
range, suggests that their effect becomes negligibly small 
at beam energies larger than ~ 500 MeV [521 . 

In conclusion, we believe that the approach described 


in this paper provide a viable and promising scheme 
for the development of a unified treatment of processes 
involving one- and two-nucleon currents, applicable in 
the kinematical region relevant to neutrino oscillation 
searches. Therefore, it may in fact be regarded as a step 
towards the new paradigm advocated by the authors of 
Refs, pi [[) . 

We are deeply indebted to A.M. Ankowski, J. Carlson, 
S. Gandolfi, C. Mariani, S. Pieper and R. Schiavilla for 
many illuminating discussions. The work of OB and NR 
was supported by INFN under grant MANYBODY. The 
work of AL was supported by the U.S. Department of 
Energy, Office of Science, Office of Nuclear Physics, un¬ 
der contracts DE-AC02- 06CH11357. Under an award of 
computer time provided by the INCITE program, this 
research used resources of the Argonne Leadership Com¬ 
puting Facility at Argonne National Laboratory, which 
is supported by the Office of Science of the U.S. Depart¬ 
ment of Energy under contract DE-AC02-06CH11357. 
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